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We consider quantum kinetics of an open quantum system in the presence of periodic fields 
designed to suppress the internal evolution and shield the system from generic low-frequency envi- 
ronment (refocusing or dynamical decoupling in application to multi-qubit systems). Assuming that 
the refocusing has order K, that is, for frozen environment the cumulant expansion of the evolution 
operator over the period r begins with the term ~ r^^^, we trace the associated cancellations in 
the kernel of the quantum kinetic equation in the Floquet formalism and characterize the remaining 
decoherence processes. 



I. INTRODUCTION 

Evolution of a quantum system subject to external 
time-dependent fields is a well studied problem that goes 
all the way to the origins of quantum mechanics. How- 
ever, driven dynamics in conventional atomic physics 
rarely involves as intricate dynamical interference pat- 
terns as those that occur, e.g., in multidimensional nu- 
clear magnetic resonance (NMR) experiments and other 
applications of coherent control (CC) where precisely 
shaped and timed signals are used to steer the quantum 
evolution of the system. One such control metho(ii*2ii 
originally developed in NMR is a pulse-based technique 
known as dynamical recoupling (also, "bang-bang" , in 
the case of hard, i5-function-like pulses). In the sim- 
plest setup, the system is a collection of individually- 
controlled weakly-coupled parts (e.g., qubits). Individual 
qubits undergo a rapid forced precession, while the over- 
all long-time evolution of the system is governed by the 
effective Hamiltonian averaged over their precession. For 
example, the interaction Jaf (jf between the two qubits is 
cancelled on average if one of them is rapidly precessing 
around the x-axis. Such cancellation of the quantum evo- 
lution is called "dynamical decoupling" or "refocusing" ; 
it is obviously related to the spin-echo experimenli^. 

An exciting thing about dynamical decoupling is its 
universality: one does not need to know the magni- 
tude of the interaction precisely to cancel it. More- 
over, with sufficiently fast pulse rate one can also can- 
cel the evolution due to slowly-varying external pertur- 
bations, in effect suppressing the decoherence caused by 
the environment^i^i'''i^i^i^°i-^^i^^i^^i-^"^i^^. General analysis 
of such methods was limited to numerics and/or the ide- 
alized (5-function-like "hard" pulses whose duration was 
either ignored or assumed small. In real experiments (es- 
pecially in solid-state systems), very short pulses are im- 
practical because they tend to couple to degrees of free- 
dom in wide spectral range which leads to signal distor- 
tions and heating—. Besides, from the previous studies it 
is hard to judge whether high-order refocusing sequences 
involving more intricate cancellation (dynamical destruc- 
tive interference) of the quantum dynamics have a real 
advantage in suppressing the decoherence. 



The goal of this work is to construct a general the- 
ory of quantum kinetics of open systems in the presence 
of periodic refocusing fields designed to suppress the in- 
ternal evolution and decouple the system from outside 
degrees of freedom. We describe the evolution of the 
density matrix of such systems, with weak internal and 
thermal bath couplings, in the approximation of mas- 
ter (quantum kinetic) equation in the Floquet formalism. 
The kinetics of the system is treated in a non-Markovian 
approximationiiiiSiiSiSS, which is essential to trace the 
decoherence suppression with the bath "slow" on the 
scale of the driven dynamics. 

Our analysis begins with the assumption that the 
period-r control fields provide ovAev-K refocusing for the 
system with the bath frozen. That is, if we replace the op- 
erators for the external degrees of freedom by c- numbers, 
the cumulant expansion of the evolution operator in pow- 
ers of thus modified internal hamiltonian begins with the 
terms of order ~ r^+^. The control fields are assumed 
to be stong and we treat them exactly. We trace the 
cancellations associated with refocusing in the kernel of 
the quantum kinetic equation (QKE) which describes the 
dissipative dynamics of the density matrix of the system 
with bath present, and characterize the remaining de- 
coherence processes. We illustrate the general analytic 
results derived for orders K < 2, which is the accuracy 
of the employed QKE, with the numerical simulations of 
dynamics of a single spin in the presence of a classical 
fluctuating magnetic field. 

Our results can be summarized as follows. Generally, 
for a weakly-coupled system weakly interacting with slow 
degrees of freedom (thermal bath) with the correlation 
time To, the decoherence is due to dissipative processes 
(resonant decay) which create excitation(s) in the envi- 
ronment, as well as reactive processes which result in de- 
phasing, or scrambling of the phase of the system. The 
associated decoherence rate is proportional to the square 
of the coupling matrix element and the correlation time 
To, see Eq. (|30|l . As a result of forced precession of the 
system caused by the control fields, the effective environ- 
ment seen by the system is modulated out of resonance, 
which may entirely suppress the state decaj»2L^. Only 
first-order [K = 1) refocusing is necessary to achieve such 
an effect. We show that, in addition, the rate of reactive 
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processes is reduced by a factor of t/tq [Eq. (03], where 
the period of the refocusing sequence, r, is assumed to be 
smaller than tq. With second-order refocusing, K = 2, 
the decoherence rate is additionally suppressed [Eq. H44I) ]. 
and with time-reversal invariant bath coupling may even 
become exponentially small in this parameter (in which 
case it will be determined by terms of higher order in 
bath coupUng, beyound the accuracy of our QKE). 

In addition to the decoherence rates which character- 
ize the exponential decay of quantum correlations with 
time, we also analyze the corresponding prefactor, which 
determines the initial decoherence^^. While for generic 
refocusing sequences with K > 1 initial decoherence is 
quadratic in r and does not scale with the thermal bath 
correlation time tq, we show that for symmetric pulse 
sequences it is reduced by an additional power of (t/tq). 

Our results extend the well-established theorySi of the 
kinetics of few-level systems in r.f. field to cases involv- 
ing a more intricate cancellation (dynamical destructive 
interference) of the quantum dynamics characteristic of 
higher-order refocusing sequences. They put in a per- 
spective the previous analyses of decoherence in the pres- 
ence of hard-pulse sequences^4S4LS42ii24iiii24iii4ii&, estab- 
lish a firm basis for future studies of decoherence scaling 
in large driven qubit systems (with and without long- 
range coupling due to long-wavelength phonons which 
may be correlated^'^ between distant qubits, contrary to 
a common assumption in the quantum error-correction 
theory--), and present an efficient alternative to opti- 
mum control techniques based directly on the master 
equationSLSSiSa. 

Some of the results regarding decoherence suppression 
in the presence of higher-order refocusing sequences were 
announced previously^^, as a justification for developing 
a technique for designing higher-order pulses and pulse 
sequences. 

II. PROBLEM SETUP 

Hamiltonian of the system. We consider iV-level 
open system with the Hamiltonian 

H = Hc{t) + Hs + Hs,, + H^, (1) 

where the oscillator bath Hamiltonian Hh = X)/^ ^m'^Ji'^a' 
has the usual form, while the control Hamiltonian 

^c^^E^"W^- (2) 

a 

the system Hamiltonian 

-f^s = 2 X/ "^oi^a, (3) 

a 

and the system-bath coupling Hamiltonian 



are expressed in terms of N x N Hermitian matrices E^, 
a — 0, ... , N'^ — 1, normalized so that Tr(I](jE^) = NSap. 
It is convenient to specify explicitly Eq = 1, choose the 
remaining matrices traceless, and define the algebra via 
the commutators and the anticommutators, 

[E„, E^] = 2iC2^E^, {E„, E^} = 2i?^^E^. (5) 

For example, for a single qubit (spin), N = 2, we can 
choose Eq, = ctq, a = 0, ... 3, in terms of the unit matrix 
(To = 1 and the three Pauli matrices, in which case the 
net coefficients B,y{t) — Vv{t) + J„ + b^, v = 1, 2, 3, can 
be interpreted as the components of the time-dependent 
magnetic field acting on the spin. 

Similarly, for two-qubit system the full set can be cho- 
sen in terms of the direct products ai^aj, i,j = 0, 1, 2, 3. 
In this case, the coefficients in front of the single-spin 
operators crii, = a^(^ ao and a2iy = cro(^ a^, can be inter- 
preted as the components of the magnetic fields acting 
on the corresponding spin, while the two-spin operators 
Ci, Up describe spin couplings. 

The same matrices will be used to parametrize the den- 
sity matrix of the system, 

P= ^(So + 5]i?aE„). (6) 

a>l 

The normalization is chosen so that Trp ~ 1. Also, for 
a pure state, = p, we have = (Ra)^ ^ N — 1 
(summation implicit), while for fully mixed state R^ = 0. 

Generally, only few of all iV^ allowed terms are ex- 
pected to be present in the Hamiltonian. Particularly, 
for n-qubit system with N = 2"^, it is common to have 
single-qubit controlling fields Viu, I — 1, ... ,n, while the 
couplings (both intrinsic ones, J^, and bath, ba) can be 
one-, two-, or multi-particle as physically appropriate. 

Dynamical decoupling in a closed system. Here 
we consider an auxiliary control problem for the sys- 
tem Q with the thermal bath operators ba in Eq. Q re- 
placed by constant numbers, which in effect modifies the 
coefficients in the system Hamiltonian (jSJ. The control 
goal is to suppress the unitary evolution with thus modi- 
fied hamiltonian Hg as fully as possible. Unless the fields 
Va{t) are specified to exactly compensate the Hamilto- 
nian Hg (which is never practical), the refocusing can 
only be achieved at some discrete set of time moments 
to — 0, ti = T, . . . . The unitary evolution over the 
refocusing interval t is commonly analyzed in terms of 
the effective Hamiltonian theory, a perturbative scheme 
based on the cumulant (Magnus) expansion for the evo- 
lution operatOB^SiSi. The expansion is done around the 
evolution in the applied controlling fields [Hamiltonian 
Hc{t)], while the system Hamiltonian Hg is treated per- 
turbatively. Obviously, this implies that the controlling 
Hamiltonian dominates the evolution. 

Explicitly, consider the evolution operator U{t), 

U{t)^-i[Hc{t)+Hg\U{t), C/(0) = 1. (7) 
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As usual, the time-dependent perturbation theory is in- 
troduced by separating out the bare evolution operator, 

Uit) = Uoit) R{t), Uo{t) = -tHcit) Uo{t). (8) 

Then, the operator R{t) obeys the equation 

m = ^tHsit)R{t), Hsit) = C/t(i) Hs Uoit), (9) 

which can be iterated to construct the standard expan- 
sion R{t) = 1 + Ri{t) + R2[t) + ... in powers of (<i?s), 



and the first-order refocusing condition, Ci(t) = 0, is 



Rn+l{t) = -tHs{t)Rn{t), i?o(t) = 1. 



(10) 



The standard Magnus expansion is readily obtained by 
integrating Eqs. H10|l formally and rewriting the result in 
terms of cumulants, 



R{t) = exp{Ci{t)+C2{t) + ...), 



dtiHsih), 



(11) 
(12) 
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C^it) = -- / dt2 / dh[Hsih),Hsit2)],---(l3) 
^ Ja Jo 

Generally, the term Ck contains a fc-fold integration 
of the commutators of the rotating-frame Hamiltonian 
HsiU) at different time moments ti and has an order 
(tHs)''- For a given Hamiltonian Hs, order- K refocusing 
is characterized by vanishing of the terms Ck of order 
k < K a,t the given time moment t = t. This is equiv- 
alent to the condition Rkir) — ior k — 1, . . . , K. The 
latter matrices can be efficiently evaluated numerically, 
which gives a systematic method for analysis and opti- 
mization of the controlled dynamics in high orders of the 
cumulant expansioni^. 

In terms of the matrices S^, the unitary transforma- 
tion generated by the control fields amounts to a rotation. 



c/o(t)i]at/t(t) = Qa^pit)^:^, 



(14) 



where the matrix Q{t) is orthogonal, Q'' = Q^^- We 
assume the matrix to be periodic with the refocusing pe- 
riod T, Q{t) = Q{t -\- t), which will be referred to as the 
"zeroth order" refocusing condition. This condition is 
non-trivial; it does not reduce to the periodicity of the 
control fields Va (t) . 

The periodicity of the real- valued matrices Q{t) implies 
the Fourier expansion with the frequencies f2„i = 2nm/T, 



Q{t) = y Arae'''"^\ A^^ = A*^. (15) 
From orthogonahty, Q{t)Q*-'{t) = i, we have 

J2^AkAll_,^S,n,oi- (16) 

With these definitions, it is easy to rewrite the first two 
refocusing conditions in algebraic form. Specifically, 

iCi(T) = ijaS^g J dtQapit), 



[Ao]i3aJa = 0, or juSt if,'' J = 0, 



(17) 



where in the second form of the expression we treated 
the coefficients Jq, as a column vector. 

Performing the double integration in Eq. (|13|l in the as- 
sumption that the first-order refocusing condition is sat- 
isfied, we have for the second-order refocusing, C2(t) — 0, 



a 



mj a/3 



0, 



(18) 



where the coefficients C^^ define the commutators, see 
Eq. ||SJ). We note that the sum in Eq. H18|l is antisym- 
metric with respect to indices a, [3, and an analogous 
condition with the symmetric coefficients [which de- 
fine anticommutators in Eq. is trivially satisfied. 



III. QUANTUM KINETICS. 

Quantum kinetic equation in rotating frame. In 

this work we consider slow (on the scale of the refocus- 
ing period r) environment, which makes it necessary to 
consider quantum dynamics of the system outside the 
commonly used Markovian approximation. We write the 
master equation as 



17.18.19.20 



m = -i[{H^{t)),p{t)\ 
■t 

dt' Trb[(5i?i(t), [(5i/i(i'),p(<')Pb]], (19) 



where (Hi(t)) = Trb (i?i (i) Pb) and 5Hi{t) = Hi{t) - 
{Hi{t)). Here Hi{t) is the interaction representation 
of the perturbation Hamiltonian Hi = Hs + Hsh [see 
Eq. in the rotating frame generated by the con- 
trol and the thermal bath parts of the Hamiltonian, 
Hq = Hc + Hh, and the bath is assumed to be in thermal 
equilibrium, pb = exp(— /3i7b)/^, Z = Trb exp(— /?_ffb)- 

With the definitions @, Q the average perturbation 
Hamiltonian in the first term of the QKE H19() is given 
just by Eq. ®, 

(i/i(<)> = Hs{t) = \[Q''{t) J]„E„, (20) 
while the corresponding fluctuating part is 

5Hi{t) - Hsb(t) - \Q'§^{t)h^{t) (21) 

where the oscillator fields in the interaction represen- 
tation ha{t) are given by Eq. with the replacement 

The second term in the r.h.s. of the QKE (|19|l is eval- 
uated in terms of two correlators, 

Tc.p{t-t') - Trb{6a(0&/3(i')Pb}, (22) 
^o.0{t-t') = Trb{6a(t)pb&/3(i')}, (23) 
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which in turn can be conveniently expressed in terms of 
the spectral coupling matrix (function) 

^"M-)-?E^^(---m)- (24) 

The fastest response time of the environment is charac- 
terized by the largest frequency of an oscillator present in 
the system. We will introduce the cut-off frequency uJc, 
such that Fai3{uj) is only non-zero for uj < ujc- In addi- 
tion, we will characterize the bath with a possibly slower 
"correlation time" tq, which describes the width of typi- 
cal features of the spectral coupling function Faf3{uj). 

Operators ba{t) are Hermitian, thus T^{t) — J-{—t). 
An explicit calculation gives !Fa(3{t) = ^afsi^)^ ^-"^^ 

/•OO J 

Ht) = J ^[F{uj){n^ + l)e'^' + F*{uj)n^e-''^']. (25) 

It is convenient to split this correlator onto real and imag- 
inary part, :F{t) = Ti{t)+T2{t). The QKE ^ becomes 

m = -'i^Q':p{t)jp[^c„p{t)\ 

-\j)it' [Q'\t)T,{t - t')Q{t')U [So, Pit')]] 
-'-j'dt' [Q''{t)p2{t ~ t')Q{t')U [So, {S^, p{t')]].l2Q) 

This can be further simplified with the help of the defi- 
nitions 10, 10): 

R-tit) = Q':^p{t)JpCl,Rs{t) 

+ / dt' [Q''{t)P,{t - t')Q{t')UCl^,cfMt') 
Jo 

+ I dt' [Q''{t)Mt - t')Q{t')UCl^,B'^pMt') 
Jo 

+ fdt'[Q''{t)Mt-t')Q{t')]apCl^, (27) 
Jo 

where the last term comes from the (time-independent) 
part of the density matrix © proportional to Sp = 1. 
This term is responsible for establishing the equilibrium 
at large t (dynamical equilibrium with refocusing). 

We note that the structure of the QKE [in particular, 
Eqs. (|25|l . (|27|l ] remains the same even if the nature of 
the thermal bath coupling is changed (e.g., by adding 



non-linear oscillator couplings) as long as the bath re- 
mains in thermal equilibrium. In such cases, the only 
change would be a renormalization of the average Hamil- 
tonian (|20|l and of the coupling matrix (|24|) . 

Kinetics in the absence of control. In the absence 
of refocusing, Q{t) — 1, the QKE does not depend 
on time explicitly, and it can be solved with the help of 
the Laplace transform [denoted with tilde, f{t) f{p)], 

pR^ip) - R^iO) = n^sip)R5ip) + ^2Mp)CIi3, (28) 
where the kernel 

= + ci^, (.Fi,o^(p)c^; + :f2Mp)b^5) (29) 

incorporates the first three terms in the r.h.s. of Eq. H27|l . 

The dissipative dynamics of the system is defined by 
the singularities of the matrix [Ip — II(p)]"^, whose lo- 
cation on the complex plane p determine the spectrum 
of the decoherence rates. For long-time dynamics, only 
the singularities close to the imaginary axis are relevant. 
With both the intrinsic interactions J and bath couplings 

weak on the scale of the bath correlation time tq, a 
good accuracy can be obtained by setting p ^ in the 
QKE kernel II(p) [Eq. (j^ ]. which is equivalent to the 
Markovian bath approximation. Then, the real parts of 
the eigenvalues of the matrix 11(0) will determine the 
spectrum of the decoherence rates. If we bunch together 
all processes causing the evolution of the density matrix, 
the maximum decoherence rate can be estimated as 

To ~max[j, A(0)ro], (30) 

where A(t) = ||.F(i)|| is a norm of the correlator ma- 
trix (|25|l . and tq is the bath correlation time introduced 
below Eq. (EH). 

Quantum kinetics in Floquet formalism. The 

full kinetic equation (|27|) in the presence of refocusing 
can also be analyzed with the help of the Laplace trans- 
formation, but in this case the structure of the solution is 
complicated by the presence of the time-dependent rota- 
tion matrices (5*''(t), Q{t'). Assuming these are periodic 
( "zeroth-order" refocusing condition), we use the expan- 
sion idll) to obtain [cf. Eq. if^ ] 

pR^{p) = '^T\,nns{p)Rs{P ~ i^ni) +ry{p), (31) 

ni 

where the kernel with the frequency transfer is 



I 

TT ^/Itn Tr^7 ,r^7 /• dt^ / [A^-^>-?^l(^)^m']a/3 ' [^^-„-?^2 (c^)^m'] a/3 ' \ , . 

m.' ^ 
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and the second term in the r.h.s. of Eq. H31|) is 



rip) 



,J 27r p 



(33) 
(34) 



The obtained expression H31|l is a set of functional equations for the Laplacc-transformcd matrix elements R-yijp) of 
the density matrix We iterate these equations to obtain a formal series in powers of n,„, 

jpi) r{p2) ^ n„j,(p) n (pi)n {P2) r{p3) , 



mi 



Pi ^ P 

mi ,m2 



P2 



{m} 



Pi 



P2 



P3 



(35) 

where p„ = p — if2m„- Generally, long-time behavior corresponds to small values of p, while that near the ends 
of the refocusing interval is governed by values of p close to any ifli, resulting in an asymptotic decomposition, 



Ra{t) — J2l 



. We will analyze the terms with different I separately, beginning with / = 0. 



To evaluate the dynamics around a given frequency fJ;, we need to carefully account for the terms singular near 
p = ifli. To this end, denote the sum of all non-singular (for small p — ifli with given /) terms connecting the terms 
with the denominators p' = p — iflii and p" = p ~ iflin , 



u\,\i„{p) = Uin_i,{p) + y nm,-i'(p) + > nm,-i'{p) 

Pi — Pi P2 



Note that this definition implies 



mi ,m2^/ 



nl''V'W = nl°L,r-,(p-*f^O- 



The entire series near p = I ^ 0, can now be written as 



[p-tni-uf\ip) 



.•(p-^^70 + E^!!L(p) 



r{p - inm) 



(36) 



(37) 



(38) 



while near p = it is 

i?[°i(p) = (p-ri[°i(p))"'ko)+E 



Sm(p) , n[,°[„(p) ^ ^o}miP) S,n'ip-i^ni) 



p - iVlm p~i^r. 



(39) 



where Eq. H33(l for r{p) was substituted for completeness. 

The analysis of the obtained expressions is dramati- 
cally simplified with at least order-one refocusing, as long 
as the couplings Jq and the bath couplings jF{t) are weak 
on the refocusing time scale r, which is assumed to be 
short on the scale of the bath correlation time tq, r <C tq. 

Indeed, the universal order-one refocusing condition, 
Eq. (|17|) , implies the disappearance of the average Hamil- 
tonian regardless of the specific values of the couplings 
Jq. Thus, in the kernel rim(p) [Eq- 1221] with m = 0, 
the first term disappears completely. Furthermore, if we 
assume the set of the fluctuating fields in the bath 
coupling hamiltonian is the same as that of the con- 
stant parameters Jq in the system hamiltonian ((JJ, the 
terms with resonant denominator (m — m' = 0) inside 



the u) integrals in Eqs. (|32(l . H34(l will also be suppressed. 
Then, for small |p|, the resonant contribution to these 
expressions will be limited by > £7, which by assump- 
tion is far out in the tail region of the spectral coupling 
function (|24|l . The remaining non-resonant contributions 
can be calculated by expanding in powers of to under the 
integrals. 

In particular, the spectrum of the dissipation rates is 
determined by the positions of the singularities of the 

QKE resolvent, (p - no°o(p))~^ [see Eq. At small 

coupling these are determined by the eigenvalues of the 
matrix n[{'Q(p = 0) [Eq. (|36l) ]. To quadratic order in pow- 
ers of perturbing Hamiltonian, with the help of Eq. I|32|) . 
we have 
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m^^O fe>0 



k+1 



mj a/3 



/35 



fc+1 



Da 



r 



(40) 



where JF^''^ (0) is the A:-th derivative of the correlator J^{t) 
[Eq. (|2SJ)] evaluated at t = 0. The corresponding maxi- 
mum decoherence rate with order-one refocusing (which 
is determined by reactive non-resonant processes) can be 
estimated as 

Ti - max[jV, A(0)t], (41) 

where A(i) = ||:^(t)|| was defined below Eq. Here 
the first expression comes from the first term in Eq. H4U|) 
and originates from the non-compensated evolution due 
to the system Hamiltonian 10), while the second term 
is an estimate of the leading order of the derivative ex- 
pansion in Eq. (|40|l . The presence of the instantaneous 
correlator can be interpreted as an effect of nearly static 
fluctuations of the coefhcients Ja due to the presence of 
the bath. Comparing with the corresponding expression 
in the absence of refocusing, we note that already with 
first-order refocusing the decoherence rate is reduced as 
long as the refocusing is fast enough, t/tq <C 1, Jt <^ 1. 

With the second-order refocusing , Eq. (|18|l . all evolu- 
tion to quadratic order in J should be compensated. To 
demonstrate this cancellation explicitly for the first term 
of Eq. (gni with p = 0, denote 

A^., ^ ^ M^^ = -M. (42) 

% \ U — it Yl 

Then the second-order refocusing condition (|18|l implies 

Ma^^a^p = Mc.0{Blp + ^Clp)l:^ = 0, (43) 

while the first term in Eq. (|40|l was obtained from the 
double commutator, 

Clearly, the first two terms in the corresponding product 
with Map are zero due to the refocusing condition (|43|l . 
while the remaining two terms cancel each other due to 
the antisymmetry of the matrix M . 

The cancellation works essentially the same way for 
the terms involving symmetric matrices, even and odd 
derivatives of the real and imaginary parts of the corre- 
lator matrix T(t) respectively, and 
(again, we use the assumption that "frozen" bath fluctu- 
ations are refocused). Thus, under most general condi- 
tions, the leading order term in the derivative expansion 
will be given by J^2(0), which gives 

r2 - A2(0)t, (44) 



where IS.2{t) = \\T2{t)\\ is defined in analogy with A(t) 
but involves only the imaginary part of the correlator 
jF{t). Formally, this term is of the same order as that 
remaining after first-order refocusing, Eq. (|41|l . We 
note, however, that this contribution represents essen- 
tially quantum effects; for temperatures not small com- 
pared with the bath cut-off scale, (3bJc ^ 1, it is expected 
to be small compared with A(i). 

In practice, the leading-order contribution to the de- 
coherence rate, Eq. H44I) . is often suppressed altogether. 
Indeed, the entire contribution of J-2{t) to Eq. (|40|l is 
identically zero for the terms involving a single spin, as 
the nested commutator-anticommutator of Pauli matri- 
ces vanishes, [cq, {(T/3, cr^}] = [the value /3 = is ex- 
cluded from the implicit summation, cf. Eqs. H2ti|) . (|27|l ]. 
For more complicated systems (e.g., involving thermal 
bath correlated across several qubits), the matrix J-{t) 
is expected to be symmetric as long as the bath is time- 
reversal invariant, that is, for real- valued spectral cou- 
pling function (|^ . F{ijj) = F*{uj). In such cases all 
terms in the derivative expansion of the second order 
contribution to decoherence rate are suppressed, which 
may result in an exponentially smaller value of r2 for 
To S> T. Such a situation where J'2{t) = and all orders 
in the derivative expansion with (0) are suppressed 
are discussed in Sec. IIVI [see Figs. El [7j. Here the sec- 
ond order contribution to the decoherence rate is seen 
to be small beyound the numerical precision already for 
To/r > 1. 

Initial decoherence. The spectrum of the decoher- 
ence rates associated with the modes around a frequency 
ill, I ^ is determined by the positions of the poles of 

the corresponding resolvent, (p — i^i — nf|(p)) , in the 
vicinity oi p — ifli. Because of the formal identity (j27J, 

^u(p) ^ n{J'Q(p — i^i), the corresponding poles are dis- 
tributed around p = if2; in an identical fashion as those 
around p = 0. As a result, at time moments commensu- 
rate with the refocusing period, t ~ t, 2t, 3t, . . . , the 
contributions with all frequencies fli add coherently, with 
the common set of decoherence rates {7} whose maxi- 
mum is determined by Eqs. 1)4111 . 1)44(1 depending on the 
order of the refocusing sequencers. 

The decoherence rates {7} determine the long-time ex- 
ponential fall-off of the refocusing accuracy. The corre- 
sponding prefactor determines the initial decoherence^^; 
it can be found as the sum of the (non-singular) ma- 
trix elements in Eqs. I)38() . ()39|l . For t sufficiently small, 
r t <C 1, the correction due to the decoherence can be ne- 
glected, and the net contribution of a sector with given 
I can be found as the sum of the residues near p = ifli. 
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For example, the total weight associated with the I = and the overall total, k = J^i '^U is 
sector can be obtained from Eq. H39|) as the coefficient in 
front of p-^R{0) at T < p < f7. 



KO = (1 - TT^) 



41(0) 



, rfn^o (p) 



dp 



p=Q 

(45) 

The weight of an ^ 7^ sector is obtained from Eq. (|38(l . 

nL"l,o(0) 



n[,l(o) 



i+E 



nL1,o (0) 



1^0 



^ <L(o) 



(47) 



To quadratic order in powers of the perturbing Hamilto- 
nian (the accuracy of the employed QKE) , and to leading 
(46) order in the derivative expansion [cf. Eq. (|40l) ]. 



frM ('n\ - Mti- n r'f M \^ I^m' jjJ*^' + •^l(0))^m-m']a/3 ' 



p — 



(48) 
(49) 



Performing the expansion to quadratic order in J and 
linear order in T and collecting various terms, we obtain 
for the overall coefficient 147|) . with the same accuracy, 

{k-i),s = [e^imuci^.c$'s. (50) 

where 

Note that this expression was derived assuming solely 
order-one refocusing, yet the constant coefficients J" give 
no contribution to quadratic order here. For a generic 
first- or second-order refocusing sequence, Q ^ q ^ 
and the initial decoherence can be estimated as 

- A(0)t2. (52) 

For sequences which produce time-reversal symmetric 
evolution, Q{t) ~ Qi—t), the Fourier components are 
real- valued, A„i = A_,„, and the sum (|51|l vanishes iden- 
tically. In such cases the initial decoherence is smaller, 
and it is determined by higher derivatives of the bath 
correlation function, e.g., 

Il^-lllsymm^ |A"(0)|t4 (53) 

for the symmetric sequence in Fig. [7| 

So far we only considered the terms ~ ^(0) which de- 
pend on the initial conditions for the density matrix. The 
remaining terms in the r.h.s. of Eqs. H39() . (|38l) provide an 
additional source of errors, as these terms are responsi- 
ble for establishing the correlations characteristic for the 
stationary state at large t; they are required to vanish at 



t ^ 0. In real-time the corresponding contributions come 
with the prefactors 1 — e*^'^'"^)*, small at commensurate 
time moments t = t, 2t, . . . because the decoherence rates 
7 are small. Additional smallness arises because the refo- 
cusing tends to average out the correlations which would 
normally appear as the equilibrium is reached. Therefore, 
we expect these contributions to be quartic, beyound the 
accuracy of the present calculation. 

IV. EXAMPLE: SINGLE-SPIN KINETICS. 

In this section we illustrate the derived general expres- 
sions on an example of a single qubit (spin) driven by 
classical fluctuating fields. Specifically, we use Gaussian 
random fields ba{t) along one (x) or all three directions, 
with the correlators 

Kit)) - 0, Kit)bfj{t')) = d^pble-'''^^'<K (54) 

where bo is the r.m.s. amplitude of the random field and 
To is the correlation time. The correlated field is gener- 
ated using the spectral filter based on fast fourier trans- 
formation (FFT) of a sequence of originally uncorrelated 
Gaussian random numbers. As a result, ba{t) are actu- 
ally periodic over the simulation interval (which is always 
long compared to tq). 

The density matrix © is described by the three- 
component vector R, i?^ = 1, whose quantum dynamics 
is described by the Bloch equation, 

R= [B(t) X R], 

where Ba{t) = Va{t) + ba{t) is the net magnetic field 
[see Eqs. Q, Q]. In terms of the vector R, the spin 
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evolution in a given classical field is a rotation; the goal 
of refocusing is to reduce the total rotation angle (f). The 
average fidelity of the refocusing, the probability for the 
qubit to remain in the original state, averaged over initial 
conditions, is equal to 1 — (1 — (cos(/)))/3. 

In Fig. we show the results of time-dependent simu- 
lations for a single spin driven by one-component random 
field with four different values of the correlation time tq. 
We plot the quantity 1 — (cos cj)) , proportional to the devi- 
ation of the average fidelity for the spin to remain in the 
same state from one, as a function of time t in units of Tp, 
a time scale equal to the interval between consecutive re- 
focusing pulses. The r.m.s. amplitude of the random field 
is the same for all curves (in fact, everywhere through- 
out this work), = 0.0355/Tp. The numerical data is 
compared with the exact analytical solution, 

rt 

1- (cos0(i)) = l-e-^'^'^*"/^, 0(i)= / b,(t)dt, (55) 
i^=6^ro^(7rV2.erf. + e--l), x^^. (56) 

For long-time asymptotics we obtain 
(cos0(i)) ^e^o^oe-^— *, 7exact=6^To(V2)'/'. (57) 




32 64 96 128 160 192 224 256 



t/Tp 

FIG. 1: "Refocusing error" (specifically, three times the de- 
viation of the average fidelity for the spin to remain in the 
initial state from one), in the absence of refocusing. The 
gaussian random field is applied along x-axis only. For differ- 
ent curves, it has the same r.m.s. amplitude bo but different 
values of the correlation time ro [see Eg. 1541 . Symbols show 
the results of simulation averaged over 900 samples of ran- 
dom field; lines show the corresponding exact results 1)55^ . 
1561 1 which for t ^ tq are also very close to the QKE result 
(not shown). See text for other notations. 

To make a connection to the quantum kinetic equation, 
we notice that in the simulations we perform the averag- 
ing over the classical random fields, instead of that over 
the quantum dynamics of the thermal bath. As a re- 
sult, the correlation matrix Tap{t — t') = {ha{t)hp{t')) is 



explicitly real- valued, J^2{t) = 0. In the absence of the 
control fields, the Laplace-transformed QKE H28I) is 

R(p) = b-n(p)]-iR(0), (58) 
l^^5{p) = Jo.e'^^-' + ^^s{p) - 5-,s^c.o.{p)- (59) 

With Ja = 0, and for the random field along the x- 
axis only, the exponent and the prefactor of the exact 
long-time asymptotics can be calculated to quadratic or- 
der in the noise amplitude by expanding the resolvent of 
Eq. H58|) around the point p — Q, 

R(p)« [p-pn'(o)-n(o)]-iR(o). 

This gives in real time [cf. Eq. H57|) ] 

(COS0)-. (1-62,2)-!^-,*^ 7=T^^ 

The simulations with refocusing were performed 
using a symmetric length-8 pulse sequence "8p" 
(XYXYYXYX), as well as a set of "concatenated" pulse 
sequencs, "4c" (XYX Y), "8 c" (XYXYXYXY), "16c" 
(XYXYXYXYXYXYXYXY), etc, where X is a 7r-pulse 
along the x-direction, X is a negative-7r pulse, and the 
longer sequences are obtained recursively by ramping 
the signs of the pulses. This concatenation procedure 
is somewhat similar but differs from that used in Ref. "sj 

We used the Gaussian pulses'^'^, as well as the first- 
and second-order self-refocusing 7r-pulses, Sl and Ql re- 
spectively, designed by the authors previously Pulses 
Sl, L = 1,2 are analogous to the first-order Hermitian 
pulses'^^ but they were constructed so that the ampli- 
tude of the signal (along with the derivatives up to 2L- 
th) turn to zero at the ends of the interval of the du- 
ration Tp. Pulses Ql, L = 1,2 are similarly designed 
one-dimensional second-order self-refocusing pulses. 

The order of the sequences with the particular pulses 
and for different directions of the applied constant field 
arc listed in Tab. HI The sequence 8p has the duration 
T = 8Tp, and so the Fourier expansion of the correspond- 
ing evolution operator Q{t) starts with the frequency 

— 2tt/{8tp). Similarly, for sequences 4c, 8c, . . . , the 
Fourier expansion starts with 27r/(4rp), 27r/(8rp), etc. 
However, due to the structure of these sequences, the 
low-frequency Fourier coefficients for sequences in the nc 
series with larger n turn out to be very small numeri- 
cally, and scaling as ^ fl^, as illustrated in Figs. [3 EJ 
^ El As a result, for relatively fast ffuctuations (small 
To) the long-time refocusing accuracy for these sequences 
can be substantially better than that of equal or shorter 
ordinary sequences (compare the slopes with tq = 2^Tp 
in Figs. El and [Tj). We also note the suppressed high- 
frequency tail of the spectra in Fig. El which illustrates 
the advantage of the pulses Q2 designedi^ specifically for 
reduced spectral widthS. 

Figs. El El show the refocusing error, 1 — (cos 0) with 
the refocusing pulses present as described in the captions. 
The amplitude of the fluctuating field 6^ (t) (along one or 
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16c 


2 


4 


6* 


2 


6* 


8* 


1 


1 


1 
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2 


4* 


6" 


4* 


8" 


> 10 


1 


1 


1 


64c 


2 
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6" 


4" 


> 10 


> 10 


1 


1 


1 


8p 


1 


1 


3 


1 


1 


3 


1 


1 
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TABLE I: Order of refocusing sequences (rows) with differ- 
ent pulse sfiapes (columns) , depending on tfie direction of the 
applied constant field. The order values listed represent the 
number of cancelled terms in the cumulant expansion (lllll 
of the evolution operator with the bath variables replaced 
by c-numbers. "G" stands for Gaussian pulses—*—, Sl and 
Ql, L = 1,2 are first- and second-order one- dimensional self- 
refocusing pulses respectively with up to 2L — 1 derivatives 
vanishing at the ends of the interval—. The superscripts 
"*" or "**" denote that the first non-vanishing cumulant is 
"small" or "very small" numerically (smaller by some two and 
four orders of magnitude respectively compared to what is ex- 
pected from naive scaling). The expansion was done numer- 
ically keeping 10 orders in the time-dependent perturbation 
theory as explained in Ref. 17^ . See text for definitions of the 
sequences. 

three directions) are the same as for data in Fig. ^ but 
the vertical scale here is reduced by some two orders of 
magnitude. This totaUy hides the curvature of the plots a 
few correlation times away from the origin, which allows 
a linear fit, 

1- {cos (j)) =^ A + Bt/ro. (60) 

The coefficients represent the initial decoherence propor- 
tional to the intercepts A with the vertical axis, and the 
decoherence rate proportional to the slopes B. We note 
that the random field used in the simulations is periodic 
with the period T = 256rp; as a result the overall error 
is almost entirely compensated towards the end of the 
simulation interval. Respectively, only the data further 
than At = 3tq from the ends of the interval was used in 
the fits. 

The pulse shapes and the fluctuating fields chosen for 
simulations in Figs. and [7| are such that the pulse se- 
quences provide at least second order refocusing. With 
solely classical correlations, ^^(t) = and J-'i{t) sym- 
metric, the decoherence rate is expected to go down dra- 
matically with increasing tq. This is exactly what is seen 
in Figs.iniEl already at tq > r the real-time graphs look 
almost horizontal and the corresponding slopes B scale 
down rapidly with increasing tq , so that they become too 
small for the numerical precision of the calculation. 

While the sequence 8p is explicitly symmetric with re- 
spect to the origin, the sequences nc are not. As a result, 
with slow fluctuations, tq/t ^ 1, the initial decoherence 
for the latter sequence tends to a constant value, as can 
be seen from the intercepts A in Fig. On the other 
hand, the intercepts tend to be much smaller in Fig. \7\ 




Sl Tp/(27t) 



FIG. 2: The Frobenius norm of the matrices Am [Fourier 
transform of the evolution matrices Q{t), see Eq. USU I with 
frequencies Qrn = 2nm/Tp for sequences nc, n = 4, 8, . . . , 64 
with Gaussian vr-pulses. The width of a pulse is chosen to 
be O.OSTp so that the discontinuity at the ends of the interval 
is numerically negligible, which results in a steep cut-off at 
high frequencies. Vertical lines mark the spectrum features of 
the parent sequence 4c, which dominate the spectrum of all 
higher-order sequences. Thin dotted lines guide the eye with 
the slope corresponding to power laws as indicated. 




£2 Tp/(Z7I) 



FIG. 3: As in Fig. I^lfor first-order self-refocusing pulses Si. 



where the symmetric sequence 8p was used. This results 
in an excellent overall refocusing accuracy. 

We have also simulated the spin dynamics under the 
4c sequence in the presence of the fluctuating magnetic 
field along the z direction, using Gaussian, Si, and Qi 
pulses which provide 0th, 1st, and 2nd order refocusing 
respectively (not shown). The decoherence rates in the 
three cases are seen as proportional to tq, independent of 
To, and vanishing rapidly with tq ^ r, as expected from 
the analytic calculations. 
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FIG. 4: As in Figs. |21 El with second-order pulses Qi. Note a 
suppression of the low-frequency part of the spectrum com- 
pared with lower-order pulses. 




FIG. 5: As in Figs. |2l for second-order pulses Q2. These 
shapes vanish at the ends of the interval along with the first 
three derivatives, which suppresses the high-frequency part of 
the spectrum. 



V. CONCLUSIONS 

In this work we discussed the kinetics of a quantum 
system subject to a pulse-based control fields of arbi- 
trary shape. We concentrated on the simplest case of 
dynamical decoupling, or refocusing, where the only goal 
is to cancel any evolution due to intrinsic or extrinsic cou- 
plings. We solved the problem in the approximation of 
a non-Markovian quantum kinetic equation, which limits 
the accuracy to quadratic order in powers of the pertur- 
bations, but considered the evolution due to the control 
fields exactly. The equations correctly represent long- 
time dissipative dynamics. The corresponding decoher- 
ence rates and the prefactor are evaluated to second order 
in powers of the small parameter, the evolution amplitude 
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FIG. 6: Refocusing error with the fluctuating magnetic field 
along the x axis as in Fig. but now in the presence of re- 
focusing sequence 8c with pulses Qi (order 2, see Tab. QJ. 
Symbols represent data averaged over 900 realizations of dis- 
order, dashes are the linear fits [Eq. II6UI1 I for data further 
than At — 3ro from the ends of the interval. Inset shows the 
fit coefficients for sequences 4c, 8c, and 16c as a function of 
the ratio of the correlation time tq. Dotted lines on the inset 
indicate the slope corresponding to the power laws as indi- 
cated. The decoherence rate (proportional to the slopes B) is 
reduced dramatically for the correlation time tq exceeding the 
duration of the sequence, t = riTp for sequence nc. Yet the 
refocusing error does not disappear altogether because of the 
initial decoherence (proportional to the intercept A) which 
does not vanish with increased noise correlation time tq for 
these non-symmetric sequences [see Eq. 11521 1. The data on 
the inset also shows that the refocusing accuracy is not im- 
proved with the longer sequences of order above 2nd, but it 
also does not worsen even for small tq/t [the low- frequency 
harmonics Am are suppressed, see Fig. The refocusing 
errors are strongly suppressed near the end of the interval be- 
cause the fluctuating field used in the calculation is periodic 
with the period T — 256rp. 



due to the perturbation over the period of the refocusing 
sequence. 

We demonstrated that higher order refocusing se- 
quences can be very effective in cancelling the decoher- 
ing effects of the couplings to slow external degrees of 
freedom. If in the absence of control the decoherence 
rate due to the bath with the characteristic correlation 
time tq is Fq [Eq. (|30|l ]. with sufficiently fast order-one 
period-r refocusing (r ^ tq) the decoherence rate can 
be reduced by a factor of ^ {t/tq) [Eq. (|41|l ]. This re- 
duction accounts for both dissipative and reactive terms 
and is dominated by the latter, as long as the driven dy- 
namics is in the spectral gap of the thermal bath. With 
second-order refocusing, the decoherence rate is further 
reduced, as it is now determined only by the quantum 
part of the bath correlator [Eq. With the bath 

coupling time-reversal invariant, additional cancellations 
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P 

FIG. 7: As in Fig. El but for the sequence 8p (pulses Qi) 
and for the magnetic field fluctuating in all three directions, 
Eq. 1122^ . Amplitude bo along each direction and other simu- 
lation parameters are as in Fig. |S| Inset shows the linear flt 
coefficients for this seqience with the magnetic field fiuctuat- 
ing in one (8p:l) and all three (8p:3) directions. The effective 
decoherence rate for the fastest fiuctuations, tq — 4rp, is big- 
ger then those in the simulations with nc sequences. However, 
with To > T = 8Tp, the decoherence rate (slope B) again goes 
down dramatically, while the initial decoherence also scales 
down quadratically with increasing ro, resulting in a superior 
refocusing accuracy. 

are possible, which may ultimately lead to the decoher- 
ence rate (in the QKE order) smaller than any power of 
the adiabaticity parameter t/tq. 

As noted on many occasions in NMR literature, sym- 
metric refocusing sequences provide for additional can- 
cellations in the evolution operator and often provide 
superior refocusing accuracy'^^. Here we show that the 
symmetry is also crucial for reducing the initial decoher- 
ence, an effective dephasing which occurs at the begin- 
ning of the refocusing sequence. While generic first- or 



higher-order control sequences result in an initial deco- 
herence proportional to the square of the amplitude of the 
fluctuating fields, ~ Tqt^/to, [Eq. (O], with symmet- 
ric sequences this leading-order contribution is cancelled, 
which produces an additional reduction by a power of the 
adiabaticity parameter t/tq. 

We illustrated these cancellations by simulations of 
a single qubit in the presence of a classical fluctuating 
magnetic field. Our simulations suggest that using non- 
symmetric refocusing sequences of order higher than two 
does little to improve the decoherence rate of the con- 
trolled system. Unlike the formulae which target the 
scaling, the simulations also illustrate the actual mag- 
nitude of the achieved reduction in decoherence. 

In this work we concentrated on the dynamics of a rel- 
atively small quantum system and ignored the scaling of 
the decoherence rates with the size of the system. For ex- 
ample, the estimate Eq. H44|l can be rewritten as an upper 
bound on the decoherence rate, in which case it contains 
an additional factor of N, the number of levels in the 
controlled system. Further studies with specific models 
of bath coupling are needed to understand in what cases 
this scaling with N can be suppressed. Present estimates 
are useful for small, few-qubit systems, or for situations 
where thermal bath does not induce long-range correla- 
tions between distant qubits. We plan to analyze the 
scaling with the system size and the range of correlations 
in the thermal bath in a future publication. Another 
planned extension of this work is to analyze the quan- 
tum kinetics of a system and ways to reduce decoherence 
during the operation of a quantum algorithm, without 
the assumptions that the control fields are periodic. 
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